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Abstract. Apparent (radiation) radius of neutron star, Roo, depends on the star gravitational mass in quite a 
different way than the standard coordinate radius in the Schwarzschild metric, R. We show that, for a broad 
set of equations of state of dense matter, Roo{M ma , K ) for the configurations with maximum allowable masses is 
very close to the absolute lower bound on Roo at fixed M, resulting from the very definition of Roo- Also, the 
value of Roo at given M, corresponding to the maximum compactness (minimum R) of neutron star consistent 
with general relativity and condition t] SO und < c, is only 0.6% higher than this absolute lower bound. Theoretical 
predictions for Roo are compared with existing observational estimates of the apparent radii of neutron stars. 
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1. Introduction 

Measuring the spectrum of radiation from neutron star 
surface (or, more precisely, atmospheres), combined with 
knowledge of distance, enables one, assuming spherical 
symmetry, to determine total luminosity, effective surface 
temperature, and neutron star radius. Recently, such stud- 
ies have been carried out for Geminga (Golden & Shearer 
1999) and RX J185635-3754 (Walter 2001); distance from 
these isolated neutron stars was obtained from measuring 
of the annual parallax. Very recently, Chandra observa- 
tions of X-ray sources in the globular clusters (whose dis- 
tance is known with relatively good precision) were pro- 
posed (and applied) to calculate the radius of a neutron 
star in quiescence (Rutledge et al. 2001). 

Neutron star are relativistic objects, and for masses 
above solar mass their radii may be only 1.5-2 times 
larger than the gravitational (Schwarzschild) radius r g = 
2GM/c 2 = 2.95 (M/M ) km. Therefore, because of a 
sizable spacetime curvature close to neutron star, one has 
to distinguish between the "true" or "coordinate" radius, 
R, which is the radial coordinate of the stellar surface 
in the Schwarzschild metric, and the "apparent radius" 
(sometimes called "radiation radius"), Roo, as determined 
by a distant observer studying radiation from neutron star 
surface. 

In the present letter we calculate dependence of Roo 
on neutron star mass for a broad set of equations of state 
of dense matter. We discuss properties of the theoretical 



Roo(M) curves, and finally we confront theoretical cal- 
culations with recent observational estimates of apparent 
radii of neutron stars. 



2. R{M), R c 
bounds 



,(M), and their theoretical lower 



The effective surface temperature, T s , at the neutron star 
surface, is related to total photon luminosity, L 7 , by 



i 7 = 4irR 2 <7 SB T* 



(1) 



where all quantities are measured by a local observer on 
neutron star surface. Spherical symmetry is assumed. A 
distant observer ( "at infinity" ) will measure "apparent 
luminosity" , "apparent effective temperature" T s °°, 
and "apparent radius" Roo, related to quantities appear- 
ing in Eq.(||) by (Thorne 1977) 



(2) 



L» = L 7 ( 1 - JL ) = 47rJ&<7SB (TrV 
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As we will show, dependence of Roo on neutron star mass 
differs considerably from R(M); the difference, which re- 
flects spacetime curvature near neutron star, increases 
with increasing M, and becomes quite large at the max- 
imum allowable mass, M max . The curves R(M) and 
Roo(M), calculated for a broad set of equations of state 
(EOS) of dense matter, are presented in Fig. 1 and Fig. 
2, respectively. One notices, that for moderately stiff and 
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Fig. 1. Neutron star radius R versus gravitational mass 
M, for seven EOS of baryonic matter, labeled by numbers 
1-7. 1: BPAL12 of Bombaci et al. (1995); 2: EoSNlHl 
of Balberg et al. (1999); 3: FPS of Pandharipande & 
Ravenhall (1989); 4: Baldo et al. (1997); 5: Douchin & 
Haensel (2000); 6: EoSNl of Balberg et al. (1999); 7: 
EoSN2 of Balberg et al. (1999). Dotted line corresponds 
to strange stars built of self-bound quark matter (SQ1, 
Haensel et al.1986). Long dashes: hybrid neutron stars 
of dense matter with a mixed baryon-quark phase, EOS 
from Table 9.1 of Glendenning (1997). Long dashes-dot 
line: EOS with first-order phase transition to a pure kaon- 
condensed matter (Kubis 2001). Doubly hatched area is 
prohibited by general relativity and corresponds to R < 
|r g . Singly hatched area is excluded by general relativity 
combined with condition w S ound < c. In the case of stars 
built of baryonic matter, configurations with maximum al- 
lowable mass is indicated by a filled circle, and in the case 
of strange stars, built of self-bound quark matter - by an 
open circle. Shaded vertical band corresponds to the range 
of precisely measured masses of binary radio pulsars. 

stiff equations of state (M max i> 1.8 M ) without a strong 
softening at highest densities, for M J> 0.5 M Q the appar- 
ent radius R^ increases with increasing M (except for a 
tiny region close to M max ), in contrast to R(M), which 
decreases in the same mass interval. 

In Fig. 1, straight lines, marking upper boundaries of 
the hatched regions of the R — M plane, result from quite 
general physical conditions imposed on the configurations 
of hydrostatic equilibrium in general relativity. The lower 
boundary results from the condition that pressure within 
an equilibrium configuration should be finite, and can 
be expressed as R(M) > |r g (Buchdahl 1959; general 
proof can be found in Weinberg 1972). This condition can 
rewritten as R(M) > # min i(M) = 3.32 (M/M Q ) km. 



Fig. 2. Apparent radius of neutron star, R^, versus grav- 
itational mass, M. Notation as in Fig. 1. Thick long-dash 
- dot straight line corresponds to ni i n i n in m Roo cit & given 
M. 

A stronger condition is obtained if we additionally re- 
quire that sound speed within the star should be sublu- 
minal: fsound = \J dP/dp < c (such a condition is neces- 
sary, but not sufficient (Olson 2000), to respect causality 
in a fluid medium). The condition w SOU nd < c implies a 
lower bound on R at a given M (Lindblom 1983, Haensel 
& Zdunik 1989, Lattimer et al. 1990, Glendenning 1992, 
Haensel et al. 1999, Glendenning 2000). In what follows, 
we will use numerical values of the absolute lower bound 
i?min2(M) = 4.17 (M/Mq) km, as calculated in Haensel 
et al. (1999); the older values obtained in (Lattimer et al. 
1990, Glendenning 1992) are slightly higher, because of 
the assumed presence of an outer envelope (crust) in neu- 
tron star models, while the value deduced from (Lindblom 
1983) is not very precise (see Haensel et al. 1999). 

A strict lower bound on R^ (M) results from the very 
definition of R^ (Lattimer & Prakash 2001). Namely, the 
definition of R^ implies 



R 



r g y/l-r g /R 



(3) 



The right-hand-side of the above equation is a function 
of x = r g /R only. It diverges to +oo at x — and at 
x = 1. At fixed M, it has a single minimum at x = 2/3. 
Therefore, minimum value of Roo(M) is i£oo,min(M) = 
7.66 (M/Mq) km (Lattimer & Prakash 2001). 

Let us notice, that this limiting R^ for an "appar- 
ently most compact" neutron star is very close to (but 
a little smaller than) that for a maximum compact- 
ness x = 0.7081 consistent with w SOU nd < c, given by 
7.71 (M/Mq) km. However, at any M the difference is 
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only 0.6%, and therefore in practice smallest J? at a fixed 
M, consistent with w SO und < c, can be considered as cor- 
responding to smallest R 00 (M), and vice versa. 

This result can be easily understood, because x — 
0.7081 is only by 0.05 higher than 2/3 corresponding to the 
minimum of Roo- Therefore, relative difference between 
-Roo.min and the value corresponding to minimum value of 
R at given M (assuming v SOU nd < c), can be estimated 
as ~ 2.25 x (0.05) 2 = 0.6%, which is consistent with our 
exact result. 

While the subluminal (w SOU nd < c) upper bound on x 
at given M is slightly larger than 2/3, the actual maximum 
values of x for various EOS, which are reached at M max 
for these EOS, are lower than 2/3. However, if we restrict 
to medium stiff and stiff EOS, with M max <; 1.8 M©, 
then £M max — 0.6, which is only 0.06 lower than 2/3. We 
may therefore expect, that for these EOS, which actually 
constitute majority of models in Figs.l, 2, (i?oo)M max will 
be only - 2.25 x (0.06) 2 ~ 0.8% larger than i? 00)min . X This 
explains, why for these EOS the points at M max are so 
close to the Roo,mm{M) line in Fig. 2. 

As shown by Lattimer & Prakash (2001), one expects 
that for any baryonic baryonic EOS, Roo > 11 km, inde- 
pendently of neutron star mass. Our Fig. 2 confirms this 
"practical lower bound" on Roo . On the contrary, there is 
no lower limit on Roo for bare strange stars, whose size can 
be as small as hundred fermis. For strange stars covered 
with a layer of normal matter, minimum radius is reached 
at M ~ 0.01 M Q , and for a maximally thick crust it is 
about 5 — 6 km (see, e.g., Glendenning 1997). 

3. Observational bounds on 

In what follows we will briefly review observational de- 
termination of apparent radius of neutron star. We will 
restrict to cases, which seem to us most promising. In 
all cases, what one determines is actually Roo/d, where 
d is neutron star distance. Therefore, independent knowl- 
edge of d is mandatory to calculate R^ from observa- 
tional data. Generally, after fitting the spectrum of pho- 
tons emitted from neutron star, one tries to get the inter- 
val i?oo.i < Roo < Roo.u, to which the value R^ belongs 
at not less than 90% confidence level. An EOS is consid- 
ered to be ruled out, if no point on its R aD (M) curve can 
satify this condition. In practice, the condition reduces to 
.Roc < J?oo,u- Generally, conclusions from the application 
of this criterion should be taken with a grain of salt, be- 
cause of the difficulty in estimating of the error in the 
photon spectrum fitting. 

3.1. Close by isolated neutron stars 

At the 95% confidence level, results of Golden & Shearer 
(1999) imply R^ < i?oo iU = 17.6 km assuming the H at- 
mosphere, and Roo < Roo,u — 16.5 km for the black body 
thermal spectrum (which turns out to be practically indis- 
tinguishable from the Fe/Si model atmosphere spectrum). 
Therefore, the value of the upper bound on the apparent 



radius of Geminga is, fortunately, not very sensitive to the 
assumed atmospheric model. 

As for this writing, the case of RX J185635-3754 
(Walter 2001, and references therein), is much less clear. 
Using atmosphere model of the photon spectrum, one de- 
duces from numbers quoted in Walter (2001) the upper 
bound i?oo iU (atm) = 18 km at the 2cr confidence level. 
However, if one uses the black body spectrum model, as- 
suming spherical symmetry, one gets an abnormally small 
upper bound, i?oo iU (BB) = 8.4 km (Walter 2001; notice 
that we use results at the 2a confidence level). Clearly, 
these results are very preliminary, and we have still to 
wait for more reliable and less model dependent determi- 
nations of Roo for this isolated neutron star. 

3.2. X-ray transients in globular clusters 

Very recently, Rutledge et al. (2001) proposed a method of 
measuring Roo of neutron stars, observed as X-ray tran- 
sients in globular clusters. As an example, they studied 
transcicnt X-ray source CXOU 132619.7-472910.8 in NGC 
5139, fitting photon spectrum with the H- atmosphere 
model. They obtained, at 90% confidence level, Roo = 
14.3 ± 2.5 km (assuming 10% uncertainty in the distance 
to NGC 5139), which results in = 16.8 km. The ad- 

vantage of the proposed method stems from the fact, that 
for neutron stars located in globular clusters both the dis- 
tance and interstellar hydrogen column density are rather 
well known. 

3.3. Other estimates of neutron star radii 

They are mostly related to accreting neutron stars, ob- 
served as X-ray bursters (see, e.g., Titarchuk 1994, Habcrl 
& Titarchuk 1995, Burderi & King 1998, Psaltis & 
Chakrabarty 1999, Li et al. 1999a, Li et al. 1999b, and 
references therein). One has to mention a strong model 
dependence of theoretical analyses, and frequent neglect 
(Burderi & King 1998, Psaltis & Chakrabarty 1999, Li et 
al. 1999a, Li et al. 1999b) of the space-time curvature ef- 
fects, which are actually crucial for the difference beteen 
R and R^. 

3.4. Theory versus observations 

As one can see in Fig. 2, the upper bounds i?oo, u on the 
apparent radius of Geminga are consistent with theoreti- 
cal predictions for Roo(M), based on considered baryonic 
EOS of dense matter, provided the neutron star mass is 
above 0.2 M . 

The value of the upper bound i?oo,u = 16.8 km, ob- 
tained by Rutledge et al. (2001) for the transient X-ray 
source CXOU 132619.7-472910.8 in NGC 5139, is con- 
sistent with neutron star R 00 (M) curves in Fig. 2, pro- 
vided M > 0.4 M©. In the case of stiffest EOS (with 
M max > 2.2 M©), the value of i?oo, u rules out high-mass 
neutron stars with M > 1.6 -f- 2.0 M©. 
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As for this writing, attempts to estimate the apparent 
radius of RX J185635-3754 (Walter 2001) are very model 
dependent. Nevertheless, we are tempted to make a fol- 
lowing comment. Had we accepted the estimate i?oo(BB) 
for RX J185635-3754, this object could be but a low-mass 
strange star. In order to produce thermal photon spec- 
trum, this low-mass strange star should have been cov- 
ered with a layer of normal matter, because a bare quark 
surface would be a too weak photon emitter (Chmaj et 
al. 1991, Usov 2001). One can only hope that more pre- 
cise measurement of the photon spectrum for this isolated 
neutron star will liberate us from such basic ambiguities. 

4. Summary 

Detection of photons emitted from the surface of isolated 
neutron stars of known distance can result in determi- 
nation of the apparent neutron star radius, i?oo- Due to 
significant space-time curvature, dependence of Roo on on 
stellar mass is quite different from that of the standard 
"coordinate" radius R. The very definition of R^ implies 
a lower bound, obtained by Lattimer & Prakash (2001), 
of 7.66 (M/Mq) km. At any M, this lower limit is very 
close to the value of Roc corresponding to the minimum 
i?, calculated under the condition of subluminal sound. 
Simultaneously, the actual values of Roo calculated at the 
maximum allowable mass are also close to this limit. For 
moderately stiff and stiff EOS with M max > 1.8 M Q , the 
actual value of R^ (M max ) is less than one percent higher 
than the absolute lower bound on R^ at this value of 
stellar mass. 

Most reliable observational estimates of i?oo , obtained 
for Geminga and transient X-ray source CXOU 132619.7- 
472910.8 in NGC 5139, lead to upper bounds on R^, 
which are consistent with all considered baryonic EOS, 
provided the mass of neutron star is above 0.2 ~ 0.4 M©. 
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